ASYMPTOTIC DENSITY OF EIGENVALUE CLUSTERS FOR THE 
PERTURBED LANDAU HAMILTONIAN 



ALEXANDER PUSHNITSKI, GEORGI RAIKOV, AND CARLOS VILLEGAS-BLAS 



Abstract. We consider the Landau Hamiltonian (i.e. the 2D Schrodinger operator 

, with constant magnetic field) perturbed by an electric potential V which decays sufR- 

I ■ ciently fast at infinity. The spectrum of the perturbed Hamiltonian consists of clusters 

, of eigenvalues which accumulate to the Landau levels. Applying a suitable version of 

■ the anti-Wick quantization, we investigate the asymptotic distribution of the eigenval- 
f— I , ues within a given cluster as the number of the cluster tends to infinity. We obtain an 
^ ■ explicit description of the asymptotic density of the eigenvalues in terms of the Radon 

, transform of the perturbation potential V. 
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! 1. Introduction and main results 

(N ■ 

1.1. Introduction. Let 

■ f d B \^ f d B " 

be the self-adjoint operator defined initially on C|^(M^), and then closed in L^(M^). 
^ I The operator Hq is the Hamiltonian of a non-relativistic spinless 2D quantum particle 

^ I subject to a constant magnetic field of strength B > 0. It is often called the Landau 

^ ■ Hamiltonian in honor of the author of the pioneering paper |23]. The spectrum of Hq 

• ^ ■ consists of the eigenvalues (called Landau levels) \q = B{2q + 1), g G Z+ = 0, 1, 2, . . .. 

^ ! The multiplicity of each of these eigenvalues is infinite, and so 

a{Ho) = aess{Ho) = |J{A J, = B{2q + 1). 

Next, let V E C(R'^; M) satisfy the estimate 

(1.1) \Vi^)\<C{^)-^, xGM^ p>l, 

where (x) := (1 + |xp)^/^. We also denote by V the operator of multiplication by V 
in L^(]R^). Consider the perturbed Landau Hamiltonian H = Hq + V. The spectrum 
of H consists of eigenvalue clusters around the Landau levels. More precisely, we have 
aess{H) = cress(-f^o) cind so all eigenvalues of if in M \ a^ssiH) have finite multiplicities 
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and can accumulate only to the Landau levels A^. Our first preliminary result says that 
the eigenvalue clusters shrink towards the Landau levels as 0{q~^^^) for g — )■ oo: 

Proposition 1.1. Assume (11 .ip .- then there exists Ci > such that for all q G Z+ one 

has 

(1.2) a{H) n [Xg -B,Xg + B]c (A^ - CiX^'^', A, + CiX-'/^). 
The proof is given in Section |H 

Remark 1.2. (i) Obviously, the above estimate 0{Xq ) for the width of the gth 
cluster can also be written as 0(g^^/^); however, as we will see, Xq provides a more 
natural scale than q. 

1/2 

(ii) Simple considerations (see Remark 3.2 in [22]) show that the estimate 0{Xq ) 
cannot be improved: the eigenvalue clusters have width > cAg with c > (unless 

= 0). This will also follow from the main result of this paper. 

(iii) Proposition 11.11 was proven in [22J for V G C(j"(]R^). The proof we give here not 
only covers the case of more general potentials V, but also is based on different 
ideas than those of [22]. 

1.2. Main result. Our purpose is to describe the asymptotic density of eigenvalues in 
the gth cluster as g — )■ oo. Let lo denote the characteristic function of the set (9 C M. 
For g G Z+ and O G i3(M), the Borel a-algebra on M, set 

^q{0) := rank 1^-1/2^^^^ (iJ). 

The measure /ig is not finite, and not even cr-finite, but if O is bounded, and its closure 
does not contain the origin, we have fJ.q{0) < oo for g sufficiently large. In particular, 
for any fixed bounded interval [a, /3] C M \ {0} we have 

(1.3) /i^([a,/3])= J2 dimKer(i7- A/) < oo 

for all sufficiently large g. Below we study the asymptotics of the counting measure iiq 
as g — )■ oo. In order to describe the limiting measure, we need to fix some notation. We 
denote by T C the circle of radius one, centered at the origin. The circle T is endowed 
with the usual Lebesgue measure normalized so that J^du = 2ti. For oj = {u;i,ijJ2) G T, 
we denote = (— a;2,wi). We set 

~ 1 f°° 

V{uj,b) = — V{bu + tu^)dt, wgT, 6 G M. 

Thus, V is (up to a factor) the Radon transform of V. In order to make our notation 
more concise, we find convenient to introduce the Banach space Xp of all potentials 
V G C(M2^M) that satisfy (HHD equipped with the norm 



(1.4) 



V||x, = sup(x)^|y(x)|. 

xeM2 
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Using this notation, by an elementary calculation one finds 

(1.5) \V{u,b)\<C,\\V\\xAby~', beR. 
Define the measure fi by 

^{0) = ^\v-\B-'0)\, OeBiR), 

where | • | stands for the Lebesgue measure (on T x M). Evidently, for any bounded 
interval [a, (3] C M \ {0} we have /i([a,/3]) < oo. Moreover, estimate (11. 5p implies that 
IJ, has a bounded support in M, and 

(1.6) [ \tfdfi{t) < oo, V£>l/(p-l). 
Our main result is: 

Theorem 1.3. Let V G C(M^) be a continuous function that satisfies (11. ip . Then, for 
any function g G C^(]R \ {0}), we have 

(1.7) lim f g^x-^dfXgiX) = [ 0{\)dfx{\). 

1^°° Jr Jr 

Remark 1.4. (i) The asymptotics (I1.7P can be more explicitly written as 

(1.8) lim \-^/^Ti0{^/\g{H - Xq)) = — [ [ g{BV{u,b))dbdu. 

q-^co 271 J J Jr 

(ii) By standard approximation arguments, the asymptotics (11. 7p can be extended to 
a wider class of continuous functions g. Further, it follows from Theorem 11.31 that 
if [a, P]CR\ {0}, and fi{{a}) = = 0, then 

lim A-i/X([«,/3]) = M[«,/3])- 

(J— >oo 

However, the assumption /x({a}) = /i({/3}) = does not automatically hold, i.e. 
in general the measure /x may have atoms. Indeed, a description of the class of 
all Radon transforms V{uj,b) of functions V G C^(M^) is well known, see e.g. 

Theorem 2.10]. According to this description, if a G C^(M) is an even real- 
valued function, then V{u,b) := a{b), 6 G M, a; G T, is a Radon transform in this 
class. Of course, if the derivative a'{b) vanishes on some open interval, then the 
corresponding measure fi has an atom. 

(iii) If F G C^(]R^), one can prove (see [22]) that the trace in the l.h.s. of (I1.8P has a 

1 /2 

complete asymptotic expansion in inverse powers of Xg , but the formulae for the 
higher order coefficients of this expansion are not known. 
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1.3. Method of proof. Let Pg be the orthogonal projection in L^(R^) onto the sub- 
space Ker (Hq — Xql). For i > 1, let Si be the Schatten-von Neumann class, with the 
norm ||-||^; the usual operator norm is denoted by ||-||. 

We first fix a natural number i and examine the asymptotics of the trace in the l.h.s. of 
(II. 8p for functions g G C^(R) such that g{X) = for small A. We have the following 
(fairly standard) technical result: 

Lemma 1.5. For any real i > l/(p — 1), the operators {H — XqYl(^x^^B,\q+B){H) and 
{PqVPq}^ belong to the trace class and 

(1.9) TtHH - X,)%,^^B,x,+B)iH)} = Tt{P,VP,Y + o{X-/~'^''), q ^ oo. 

The proof of this lemma is given in Subsection 14.31 This lemma essentially reduces the 
question to the study of the asymptotics of traces of {PqVPqY- Our main technical 
result is the following statement: 

Theorem 1.6. Let V satisfy (II. ip and let i?o > 0. 

(i) For some C = C{Bq), one has 

(1.10) sup sup X\'^B-^\\PqVPq\\ < C\\V\\x,. 

q>0 B>Bo 

(ii) For any real i > l/(p — 1), we have PqVPq G Si, and for some C = C{Bo,i), the 
estimate 

(1.11) sup sup X^^'-''^/^^'^B-'\\PqVPq\\i < C\\V\\x, 

q>0 B>Bo 

holds true. 

(iii) For any integer £ > 1 / [p — 1) , we have 

(1.12) \im Xi'-^^^^TiiPqVPqY = — [ [ V{io,bYdbduj. 

Although in our main Theorem 11.31 the strength B of the magnetic field is assumed to 
be fixed, we make the dependence on B explicit in the estimates ( ll.lOp . (II. lip , as these 
estimates are of an independent interest (see e.g. [10]), and can be used in the study of 
other asymptotic regimes. 

The proof of Theorem 11.61 consists of two steps. In Section [2] we establish the unitary 
equivalence of the Berezin-Toeplitz operator PgVPg to a certain generalized anti-Wick 
pseudodifferential operator (\E'DO) whose symbol Vb is defined explicitly below in ( I2.27p . 
Further, in Section 12] we study this \E'DO, prove appropriate estimates, and analyze its 
asymptotic behavior as g — )■ oo. 

A combination of (II. 9p and (I1.12p essentially yields ( II. 7p for a function g G (R) such 
that ^(A) = A^ for small A. After this, the main result follows by an application of the 
Weierstrass' approximation theorem; this argument is given in Subsection 14.41 



Remark 1.7. In |22] the limit ( I1.12p was computed for i = 1,2, but the result was 
written in a form not suggestive of the general formula. 
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1.4. Semiclassical interpretation. Consider the classical Hamiltonian function 

(1.13) H{i,^) = {i + \Byf + {r,-\Bx)\ ^ := (e, r/) G M^, x := (a;, y) G M^, 

in the phase space T*M? = with the standard symplectic form. The projections onto 
the configuration space of the orbits of the Hamiltonian flow of are circles of radius 
\/E/B, where E > is the value of the energy corresponding to the orbit. The classical 
particles move around these circles with period Tb = n/B. The set of these orbits 
can be parameterized by the energy E > and the center c G of a circle. Let us 
denote the path in the configuration space corresponding to such an orbit by 7(0, E, t), 
t G [0, Tb), and set 

(1.14) {V)ic,E) = ^ r\i^{c,E,t))dt, TB = n/B. 

For an energy E > consider the set Me of all orbits with this energy. The set Me is 
a smooth manifold with coordinates c G M^. It can be considered as the quotient of the 
constant energy surface 

Ss = {(tx)GM^|H(^,x) = E} 

with respect to the flow of %. Restricting the standard Lebesgue measure of to S^; 
and then taking the quotient, we obtain the measure Bdcidc2 on Me- An elementary 
calculation shows that the r.h.s. of f ll.Sp can be rewritten as 

(1.15) — [ [ Q{BV{u,h))dhdLo= — lim -\= [ g{^/E {V){c, E)) B dcidc2. 
27r J J 27r e^oo Jr2 

The basis of this calculation is the fact that as — i- 00, the radius \/E/B of the circles 
representing the classical orbits tends to infinity. Thus, the classical orbits approximate 
straight lines on any compact domain of the configuration space. 
Given (ll.lSp . we can rewrite the main result as 

(1.16) lim ^Tr^(yA,(iJ- A,)) = ^ lim ^ / gi^E {V){c, E)) B dci dc2. 

This agrees with the semiclassical intuition. Formula f ll.l6p corresponds to the well 
known "averaging principle" for systems close to integrable ones. This principle states 
that a good approximation is obtained if one replaces the original perturbation by the 
one which results by averaging the original perturbation along the orbits of the free 
dynamics. This method is very old; quoting V. Arnold [21, Section 52]: "In studying 
the perturbations of planets on one another. Gauss proposed to distribute the mass of 
each planet around its orbit proportionally to time and to replace the attraction of each 
planet by the attraction of the ring so obtained" . 



1.5. Related results. 



6 



A. PUSHNITSKI, G. RAIKOV, AND C. VILLEGAS-BLAS 



1.5.1. Asymptotics for eigenvalue clusters for manifolds with closed geodesies. In spec- 
tral theory, results of this type originate from the classical work by A. Weinstein |37] 
(see also [D]). Weinstein considered the operator —A_m + V, where A_m is the Laplace- 
Beltrami operator on a compact Riemannian manifold A4 with periodic bicharacteristic 
flow (e.g. a sphere), and V G C(A^;M). In this case, all eigenvalues of A_m have finite 
multiplicities which however grow with the eigenvalue number. Adding the perturba- 
tion V creates clusters of eigenvalues. Weinstein proved that the asymptotic density 
of eigenvalues in these clusters can be described by the density function obtained by 
averaging V along the closed geodesies on J\4. Let us illustrate these results with the 
case A4 = S^. It is well known that the eigenvalues of — A§2 are Xg = q{q + 1), g G Z^, 
and their multiplicities are dg = 2q + 1. For V G C(§^; M) set 



where Cuj{s) G is the great circle orthogonal to u, and s is the arc length on this 
circle. Then for each g G C^(§^; M) we have 



where dS is the normalized Lebesgue measure on S^. Since can be identified with its 
set of oriented geodesies Q, the r.h.s. of (I1.17P can be interpreted as an integral with 
respect to the S'0(3)-invariant normalized measure on Q. This result admits extensions 
to the case = §" with n > 2, and, more generally, to the case where is a compact 
symmetric manifold of rank 1 (see [STJ |9]). 

In more recent works [331 EBl |35], the relation between the quantum Hamiltonian of 
the hydrogen atom and the Laplace-Beltrami operator on the unit sphere is exploited, 
and the asymptotic distribution within the eigenvalue clusters of the perturbed Hamil- 
tonian hydrogen atom is investigated. The asymptotic density of eigenvalues in these 
clusters was described in terms of the perturbation averaged along the trajectories of 
the unperturbed dynamics (i.e. the solutions to the Kepler problem). 
Among the main technical tools used in [331 EH ES] which originate from [151 131] , are 
the Bargmann-type representations of the particular quantum Hamiltonians considered, 
implemented via the so-called Segal-Bargmann transforms. In our analysis generalized 
coherent states and associated anti-Wick ^'DOs closely related to the Bargmann repre- 
sentation and the Segal-Bargmann transform appear again in a natural way (see Sec- 
tion |2|), although their role is different from the one of their counterparts in [331 l36l [35]. 
Although this paper is inspired by [371 ESI ESI ES], much of our construction (see Sec- 
tion Hj) is based on the analysis of [22]. In [22] it was proven that for V G C^(M^) the 
trace in the l.h.s. of fll.l2p has complete asymptotic expansions in inverse powers of 

1/2 

Xq . However, the coefficients of this expansion have not been computed explicitly; see 
Remark 11.71 above. 




(1.17) 




g{V{uj))dS{uj) 



EIGENVALUE CLUSTERS 



7 



1.5.2. Strong magnetic field asymptotics. It is useful to compare our main result with 
the asymptotics as i? — )■ oo of the eigenvalues of H. It has been found in [26j (see also 
[n]) that 

(1.18) Mm B-^ Ti g{H - \g) = ^ [ g{V{^))d^= [ g{t)dm{t) 

where q G C^{^ \ {0}), q e 1+, V e ^^(R^), p > 1, and m(0) := ^\V~^{0)1 
O G i3(]R). Similarly to Theorem ll.3[ the proof of fll.l8p is based on an analogue of 
Theorem 11.61 (i) - (ii) (see Lemma [2.121 below), and the asymptotic relations 

(1.19) \im B-^i:i{PqVPqY = \im B-^i:iPqV^Pq = ^ ! V(x)^rfx, g G Z+, 

with V G Cq^(M2) and i e N, close in spirit to ffT:T2D . Since Bm{[a,P]) = 
^ \Vj^\[a, [a,P] C M \ {0}, where Vb (see ([52ZD) is the symbol of the gener- 
alized anti-Wick \E'DO to which PqVPq is unitarily equivalent, we find that f ll.l8p is 
again a result of semiclassical nature. However, f ll.lOp implies that in the strong mag- 
netic field regime the main asymptotic terms of TT{PgVPgY and TrPqV^Pq coincide, 
and hence in the first approximation the commutators [V, Pq] are negligible, while fll.l2p 
shows that obviously this is not the case in the high energy regime considered in the 
present article. Hence, Theorem 11.31 retains "more quantum flavor" than f ll.l8p . and 
hence its proof is technically much more involved. 

1.5.3. The spectral density of the scattering matrix for high energies. In the recent work 
[8] inspired by this paper, D. Bulger and A. Pushnitski considered the scattering matrix 
S{X), A > 0, for the operator pair (—A + V, — A) where A is the standard Laplacian 
acting in L^(M'^), d > 2, and V G C(]R'^;M) is an electric potential which satisfies an 
estimate analogous to ( II. ip . Although the methods applied in [B] are different from 
ours, it turned out that the asymptotics as A — )■ oo of the eigenvalue clusters for S{X) 
are written in terms of the X-ray transform of V in a manner similar to fll.7p . 

2. Unitary equivalence of Berezin-Toeplitz operators and generalized 

anti-Wick \EfDOs 

2.1. Outline of the section. From methodological point of view, this section plays a 
central role in the proof of Theorem II. 3[ Its principal goal is to establish the unitary 
equivalence between the Berezin-Toeplitz operators PqVPq, q G Z+, and some general- 
ized anti-Wick ^DOs Op7(Vb) whose symbol Vb is defined explicitly in (jMB- This 
equivalence is proved in Theorem 12.111 below. The \I'DOs Opq^ introduced in Subsec- 
tion 12. 3[ are quite similar to the classical anti-Wick operators Opg"' (see [U Chapter V, 
Section 2], [2H1 Section 24]); the only difference is that the quantization Opg"' is related 
to coherent states built on the first eigenfunction (fQ of the harmonic oscillator fl2.4p . 
while Opq^, g G N, is related to coherent states built on its gth eigenfunction ipq. 
In our further analysis of the operator Op^^IVb) performed in Section [3l we also heavily 
use the properties of the Weyl symbol of this operator. Thus, in Subsections 12.11 - 12.21 
we introduce the Weyl quantization Op'", and in Subsection 12.41 we briefly discuss its 
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relation to Op^"'. In particular, we show that Opg"'(s) = Op"'(s * \E'g) where s is a 
symbol from an appropriate class, and 27r\E'q is the Wigner function associated with yj^, 
defined explicitly in f l2.13p . Therefore, the Berezin-Toeplitz operator PqVPg, q G Z+, 
with domain PqL^(M^), is unitarily equivalent to Op"'(Vb * \E'g) (see Corollary 12. ISp . 

2.2. Weyl ^DOs. Let d > 1. Denote by 5(M'^) the Schwartz class, and by S'{R'^) its 
dual class. 

Proposition 2.1. [201 Lemma 18.6.1] Let s E ^'(R^'^). Assume that s G L1(M2'^) where 
s is the Fourier transform of s, introduced explicitly in (13.11) below. Then the operator 
Op'"(s) defined initially as a mapping from the S(R'^) into S'{M.'^) by 

(2.1) (Op'"(s)m) (x) = (27r)-^ / / s f^^,e) e'^"-"'>«M(x')c/x'c/e, x G R'', 

Jr^ Jr^ V 2 / 

extends uniquely to an operator bounded in L^(R'^). Moreover, 

(2.2) ||Op-(s)||<(27r)-is|Ui(K2.). 

Some of the arguments of our proofs require estimates which are more sophisticated 
than (Q. Let r(R2'^), d>l, denote the set of functions s : R^'^ ^ C such that 

||s||r(jR2d) := sup sup \d"d^s{x,^)\ < oo. 

{a,/3eZ^ I |a|,|/3|<[|]+l} (x,0eIR2d 

Proposition 2.2. |6l Corollary 2.5 (i)] There exists a constant cq such that for any 
s G r(R2'^), d>l, we have 

||Op"'(s)|| < co||s||r(R2d). 

Further, if s G L^{M?'^), then, obviously, the operator Op"'(s) belongs to the Hilbert- 
Schmidt class, and 

(2.3) l|Op-(.)||^ = -^ / \six,0\'dxdC 

(2vr)^ J^2d 

Next, we describe the well known metaplectic unitary equivalence of Weyl \E'DOs whose 
symbols are mapped into each other by a linear symplectic change of the variables. 

Proposition 2.3. [11, Chapter 7, Theorem A.2] Let k : R^"' R^'^, d>l, be a linear 
symplectic transformation, Si G r(R^'^), and S2 := Si o k,. Then there exists a unitary 
operator U : L'^{R'^) ^^(R^) such that 

Op"(s2) = U*Op'"{si)U. 

Remark 2.4. (i) The operator U is called the metaplectic operator corresponding to 
the linear symplectic transformation k. There exists a one-to-one correspondence 
between metaplectic operators and linear symplectic transformations, apart from 
a constant factor of modulus 1 (see e.g. [20l Theorem 18.5.9]). Moreover, every 
linear symplectic transformation k is a composition of a finite number of elementary 
linear symplectic maps (see e.g. [20| Lemma 18.5.8]), and for each elementary linear 
symplectic map there exists an explicit simple metaplectic operator (see e.g. the 
proof of ^ Theorem 18.5.9]). 
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(ii) Proposition 12.31 extends to a large class of not necessarily bounded operators. In 
particular, it holds for Weyl \E'DOs with quadratic symbols. 

2.3. Generalized anti-Wick \l'DOs. In this subsection we introduce generalized anti- 
Wick \E'DOs . These operators are a special case of the \I'DOs with contravariant symbols 
whose theory has been developed in [3]. Introduce the harmonic oscillator 

(2.4) h:=-^ + x', 

self-adjoint in L^(M). It is well known that the spectrum of h is purely discrete and 
simple, and consists of the eigenvalues 2g + 1, g G Z_|_, while its associated real- valued 
eigenf unctions ipq, normalized in L^(M), could be written as 

(2-5) ^= (^y^^' 

where 

(2.6) H,(x) := (-l)'?e"'/' ~ ^) ^ ^ ^ ^' ^ ^ 

are the Hermite polynomials. Introduce the generalized coherent states (see e.g. [28] ) 

(2.7) ip,,^^^{y) ■= e'«V,(z/ -x), y G M, (x, G R\ 
so that ifq = (Pq-ofl. Note that if / G i^^(IR), then 

(2.8) (2vr)-i / \{f,^q,^,^)\'dxd^ 

where (■, ■) denotes the scalar product in L^(M). Introduce the orthogonal projection 

(2.9) Pq.,,^:=\^q,,,^){^q,,,^\:L\R)^L\R), g G Z+, {x,0^^'- 
Let s G L\R^) + L°°(M2). Define 

Op7(s) := (27r)-i / s{x,0p,;.4dxd^ 
as the operator generated in L^(M) by the bounded sesquilinear form 

(2.10) FqM,9]-=i2n)-' [ s{x,0U\V>,;^,d{9,fg;cc,ddxd^, f,geL\R). 

JR2 

We will call Opq^{s) operator with anti- Wick symbol of order q equal to s. We introduce 
these operators only in the case of dimension d = 1 since this is sufficient for our 
purposes; of course, their definition extends easily to any dimension d > 1. 
Note that the quantization Op^"' with g = coincides with the standard anti- Wick 
one (see e.g. [29j Section 24]). In the following lemma we summarize some elementary 
basic properties of generalized anti-Wick \E'DOs which follow immediately from the 
corresponding properties of general contravariant \E'DOs . 
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Lemma 2.5. |29l Section 24] H Section 5.3] (i) Let s e L°°(R2). Then we have 

(2.11) l|Opr(s)ll < MLo^m- 

(ii) Let s G L^(]R^) with £ G [1, oo). Then we have 

(2.12) I|0p7(.)||^<(27r)-i.||l.(^.). 

2.4. Relation between generalized anti-Wick and Weyl ^DOs. For g G Z+ set 

(2.13) ^,ix, = ^L,(2(x2 + e'))e-(^'+«'\ (x, G R^ 

71 

where 

(2.14) M*)-^e'J^ = ^(yLL, i.E. 
are the Laguerre polynomials. 

Lemma 2.6. For a /ixec? (x,^) G tt;e /iawe Pg-x,s, = Op'" (?"g;x',5) where Pq-x,^ is the 
orthogonal projection defined in fl2.9p . and 

(2.15) W(a^''0 :=27r^,(x'-x,r-0, {x-O^^'- 

Proof. Using the well-known relation between the Schwartz kernel of a linear operator 
and its Weyl symbol (see e.g [20l Eq. (18.5.4)"]), we find that the Weyl symbol <iq;x,i of 
the projection Pq^x,^ satisfies 



(2.16) <^^.^.^(x',e') = / e-'<Vq,xA^' + v/2)ipq,xd^' -v/2)dv. 
By dlSD and dZ?]), 

/ e-'""^' ifiq,x,d^' + v/2)ipq.x,d^' - v/2)dv = 
Jr 

(2.17) -ri^ I e™(«-«')Hg(x' + it;-a;)H,(x'--t;-x)e-("'+^"~")'/V("'-^''-")'/2dv. 

Changing the variable of integration v = 2{t + i{^ — ^')) , and bearing of mind the parity 
of the Hermite polynomial Hg, we get 

f e^^«'«')H,(x' + ^v- x)Rq{x' -^v- x)e-(^'+^'^-^)'/2g-(x'-|.-x)V2^^ ^ 
Jr 2 2 

(2.18) 2(-l)'?e-("'-")'~(«'-«)' / e~''}lq{t-{x-x'-i{^-0))^<i{t + x-x' + t{^-0)dt. 

Jr 

Employing the relation between the Laguerre polynomials and the integrals of Hermite 
polynomials (see e.g. [HI Eq. 7.377]), we obtain 

/ e-'^Rqit -{x-x'- z(e - e')))Hg(t + x-x' + t{^- C))dt = 
Jr 
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(2.19) v^2^gM2((x'-x)2 + (e-n'))- 

Putting together ( l216l) - (^A9\\ . we obtain (I2A51) . □ 

Remark 2.7. Let G L^(M) and = 1. Then the Weyl symbol of the rank-one 

orthogonal projection is called the Wigner function associated with ip (see e.g. 

{24^ Definition 2.2]). Thus, Lemma [2.61 tells us, in particular, that 27r\l'g is the Wigner 
function associated with ipg. 

Lemma 12.61 immediately entails the following 

Corollary 2.8. Let s G L^W^) + L°°(M2). Then we have 

(2.20) Op7(s) = Op"'(^g*s). 

2.5. Metaplectic mapping of the operators Hq, Pg and V. For x = {x, y) E 

M2, ^ = (^,r/)GM2, set 



(2.21) xb(x, ^) := ( ^(x - r^), -^(e - y), ^(e + y), -^{v + x) 



Evidently, the transformation xb is linear and symplectic. Define the unitary operator 
Ub : L\R^) L\R^) by 



(2.22) (^UBu){x,y):=^ / e^'f^^^^'^-'^'^yMx' ,y')dxdy' 



where 



(j)B{x, y; x', y) := B— + B^^^{xy' - yx') - x'y' 



Writing hb as a product of elementary linear symplectic transformations (see e.g. [20| 
Lemma 18.5.8]), we can easily check that Ub is a metaplectic operator corresponding to 
xg. Note that 

n o xb(x, ^) = B{e + x^), X = {x, y) G M^ $, = (^, r/) G M^ 

where "H is the Weyl symbol of the operator Hq defined in f ll.lSp . On the other hand, 
-B(^^ + x"^) is the Weyl symbol of the operator B{h ® ly) self-adjoint in L^(M^ ,^) where 
h is the harmonic oscillator f l2.4p . acting in //^(M^.), and ly is the identity operator in 
L^(Mj,). Denote by Pg = \({)q){ipq\ = pg;o,o the orthogonal projection onto Ker (/;, — 2g — 1), 
q G Z+. Applying Proposition 12.31 with k = x^, and bearing in mind Remark 12.41 (ii), 
we obtain the following 

Corollary 2.9. (i) We have 

(2.23) U*bHqUb = B{h®Iy), 

(2.24) U^PjiAB =Pg^Iy, qe Z+. 
(zz) IfV e T{R^), then 

(2.25) U*bVUb = Op^(Vb) 
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where 

(2.26) Vb(x, y; ^, v) ■= ViB-'/'ix - v), B'^/\i - y)\ (x, y- ^, v) e 

Remark 2.10. Various versions of the symplectic transformation in f l2.2ip and the 
corresponding metaplectic operator Ub in (12.22^ have been used in the spectral theory 
of the perturbations of the Landau Hamiltonian (see e.g. [IB])- Of course, the close re- 
lation between the Landau Hamiltonian Hq and the harmonic oscillator h is well-known 
since the seminal work [23] where the basic spectral properties of Hq were first described. 



2.6. Unitary equivalence oi PqVPq and Op7(Vb). Set 

(2.27) Vsix, y) = V{-B-^l^y, -B~^'^x), (x, y) G 
Theorem 2.11. For any V G /.^(M^) + L°°(m2) ^ ^ /^^^g 

(2.28) U*sPgVPqUB =p,® Op7(\/b). 

For the proof of Theorem 12.111 we need some well known estimates for Berezin-Toeplitz 
operators: 

Lemma 2.12. [25| Lemma 5.1], [TS] Lemma 5.1] Let V G L\^^), i G [l,oo). Then 
PqVPq G Se{L\R^)), and \\PgVPq\\', < f 111^111.(^2), q e Z+. Moreover, zfV e L\R^), 
then 

(2.29) TiPqVPq = ^! \/(x)cix, gGZ+. 

Proof of TheoremlKm Assume at first V G C^(M2). Then, by ([2211), (1225]), and 

(2.30) U^PqVPqUB = {Pq ® /, ) Op"' ( ) (p, ® /,). 

Let u G 5(M2). Set 

:= / u{x,y)Lpq{x)dx. 



Then we have 

{U^PgVPqUBU,u)L2(R2) = {Op'"(VB)iVq ® Uq), {ifq ® Uq)) 

^ ^ VbUvi + 2/2)/2 - r/ - (xi + X2)/2) e^[(-i-2)?+(s/i-2/2)r,] ^ 



(27r)^ 



V5g(a;i)^^q(2/i) ^qix2)uq{y2) dxidx2 dyidy2 d^drj = 
77^2 I ^B{{yi + l/2)/2 - y', - V) e^^^-^^)" X 



1 

2^ 



V5q(V + t;/2)v95(V - t;/2)e*''^'(it;j Uq{yi)Uq{y2)dr]'dr]dy'dyidy2 = 
^q{y\ 7]')VB{{yi + y2)/2 -y\r]- vy'^''''''^''uq{y,)^^dy,dy2dy'diidri' 



EIGENVALUE CLUSTERS 



13 



(2.31) 

To obtain the first identity, we have utihzed Corollary 12.91 To establish the second 
identity, we have used (12. ip . (I2.26p . and (I2.27p . To get the third identity, we have 
changed the variables xi = 1]' + v/2, X2 = r]' — v/2, ^ = y'. To obtain the fourth identity, 
we have used (12.150 - (I2.16P with ^' = y', x' = r]', and x = 0, ^ = 0, taking into account 
that "^qirj', —y') = "^qiy', i]'). To deduce the fifth identity, we have applied (12. ip . bearing 
in mind the symmetry of the convolution * Ve = * ^ qi and for the sixth identity, 
we have applied (I2.20p with s = Vb- Finally, the last identity is obvious. Now, (I2.3ip 
entails flT^ in the case V e C^{^^). 

Further, let V G L^(M^), and pick a sequence {Kn} of functions Kn G C[j"(M^) such that 
Vm -^V'm L^{E?) as m ^ oo. Then by Lemma [2.121 and the unitarity of Ub-, we have 

lim WWBPqVmPqUB-WBPqVPqUBh = 0. 

m— >oo 

Similarly, it follows from (KT2^ with £ = 1 and (KWi that 

lim \\pq ® Op7(K^,b) -pq® Opr(VB)||i = 0. 

Hence, ^TM is valid for V e L^iR"^). 

Finally, let now V = Vi + V2 with Vi e L^{R'^) and V2 G L°°{R'^). Denote by xr 
the characteristic function of a disk of radius R > centered at the origin. Then 
Vi + XrV2 e L1(R2). Evidently, 

W-limW* P,(l^l + XRV2)PqKB = UlPqVPqUB. 

R—^oo 

while fCTIl) entails 

w-limp, ® Opr ((V^i + XrV2)b) =Pq<^ Opr (Vb), 
R-^00 

which yields (I2.28P in the general case. □ 
Combining Theorem 12.111 and Corollary 12. 8[ we obtain the following 
Corollary 2.13. Let V G L^{R^) + L°°{R^) and g G Z+. Then we have 

(2.32) U*BPqVPqUB =Pq® O^^ {Vb * ^,). 

Remark 2.14. To the authors' best knowledge, the unitary equivalence between the 
Toeplitz operators PqVPq, q G Z4., and \E'DO with generalized anti-Wick symbols in 
the context of the spectral theory of perturbations of the Landau Hamiltonian, was first 
shown in [25]. Related heuristic arguments can be found in [28l[7]. In the case g = this 
equivalence is closely related to the Segal-Bargmann transform in appropriate holomor- 
phic spaces which, in one form or another, plays an important role in the semiclassical 
analysis performed in [5U 1321 ES]- Let us comment in more detail on this relation. 
The Hilbert space PqL'^{E?) coincides with the classical Bargmann space 

|/GL^(M^)|/(x) = e-^W^/^^(x), |^ + ^| = 0, x=(x,y)GM^}. 
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Then the Segal-Bargmann transform Tq : L^(]R) — )■ PoL'^{M.'^) is a unitary operator with 
integral kernel 

%:=^( -] ''^ e-^«^+*^+=^*)'-2*'+l-l')/^ X = (x, y) G M^, tE M, 

(see [Ml Lemma 3.1]). Fix q e Z+. Denote by Mq : ^^(M) -> {p^ ® 4)L2(R2) the 
unitary operator which maps u G L2(R) into B^l^^q{x)u{B^'^y), {x,y) G M^^ and by 
7?. : L^(R^) — > L^(R^) the unitary operator generated by the rotation by angle 7r/2, i.e. 
{7lu){x,y) = u{y, —x), {x,y) G R^; note that [}Z,Pq] = 0. Then we have 

To = TZUbMq. 

From this point of view the operators Tg := TUAsMq, g G N, could be called generalized 
Segal-Bargmann transforms. 

3. Analysis of Op'"(Vb * ^g) and proof of Theorem 1.6 

3.1. Reduction of Op'"(Vb * ^g) to Op'"(Vb * 5^/2^)- In the sequel we will use the 
following notations. For A; > 0, let 5k be the ^-function in R^ supported on the circle of 
radius k centered at the origin. More precisely, the distribution 5k G iS'(R^) is defined 
by 

1 /■^'^ 

:= _ / i^{k cos e,k sin e)de, G ^(R^). 
27r Jo 

Next, we denote by / the Fourier transform of the distribution / G iS'(R'^), unitary in 
L2(R<^), i.e. 

(3.1) /(O := (27r)-'^/2 / e-'^<f{x)dx, ^ G R^ 
for / G 5(R'^). 

Lemma 3.1. Let V e ^^(R^) and Bq > 0. Then for some constant C = C{Bq) one 
has 

(3.2) sup sup XI/'B-'\\Op-{Vb * v]>,) - Op-{VB * 5^^)11 

g>0 B>Bo 

<c [ (|c|^/2 + |cr)|i>(C)MC, 

(3.3) sup sup Af i?-i||Op"(VB * ^g) - Op'^iVn * 5^^)112 

g>0 -B>Bo 

<c([ {\C\' + \cn\v{0\'dC 
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The intuition behind this lemma is the convergence of to 5^/2g+T ^in appropriate 
sense as g — )■ oo. We also note the similarity between the definition of Vb * 6^ and the 
"classical" formula f ll.l4p . For brevity, we introduce the short-hand notations 

(3.4) Sq := Vb * ^'g, qE Z+, tk := Vb * 4, ke (0, oo). 

Proof. First we represent the symbols Sq, tk in a form convenient for our purposes. For 
Sq we have 



(3.5) Sq{z 



[ e'^^^qiOVBiOdC, ze 



In the Appendix we will prove the formula 

(3.6) %{C) = i-iy^qi2-'C)/2, qez+, C e 

By ([33]), dMD, and the definition f l27[3|l of ^g, 

1 



3^^^L,(|Cr/2)e-l^l^/Vs(C)rfC- 



For tk we can write 



(3.7) t,{z) = ^ [ e'^^MmWBiOdC, z G M^ 

since the integral representation for the Bessel function Jq can be written as 

(3.8) Jo(fc|C|) = 27r4(C), C e 

Thus (13. 2p (resp. (13.30 ) reduces to estimating the operator norm (resp. the Hilbert- 

Schmidt norm) of the operator with the Weyl symbol 

(3.9) 

^/ e^<(L,(|C|V2)e-l^l'/^- Jo(v/2^|C|))i>i,(C)rfC, ? G 



Sg(2;) - 1^2^(2;) 



In what follows the estimate 



(3.10) Lg(x)e-^/2- Jo(V(4g + 2)x) < ^(g-^/V/^ + g^V), g G N, x > 0, 

plays a key role. This estimate is probably well known to experts, but since we could 
not find it explicitly in the literature, we include its proof in the Appendix. 
Let us prove the estimate (13.21) . Using the estimates (12. 2p and (I3.10p . we obtain: 

(3.11) WOv'^iyB * vj/g) - Ov'^iyB * 5^)\\ < i27r)-'\\sq - t^lUi 



(27r)-^ / |L,(|C|V2)e-l^l^/^- Jo(v/2^|C|)||1^b(C)MC 

<C [ {q-'/'\C\'^' + q-'m\VBiO\dC. 

JR2 
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Recalling the definition 1^22^ of Vb, we obtain Vb{C) = and so the l.h.s. 

of (13.111) can be estimated by 

CBq-'/' [ \Cf^\V,{B^/\)\dC + CBq-' f mV,{B''^C)\dC 

= cB~'^/'q-^/' [ \cf'mc)\dc+cB-\-' f \c\'mc)\dc 

This yields flO) . 

Next, let us prove the estimate (13.31) . By (12.31) and the unitarity of the Fourier transform, 



= (27r)-i I |L,(|Cr/2)e-l^l'/'- Jo(v/2^|C|)r|V^^(C)rrfC 



Now using the estimate (l3.1Up again, we obtain (13. 3p in a similar way to the previous 
step of the proof. □ 

3.2. Norm estimate of 0^'^{Vb * 4)- 

Lemma 3.2. Let V{^) = {x)-p, p > 1, and Bq > 0. Then 

sup sup kB-'^/^\\0p'^{VB*6k)\\ < oo. 

fc>0 B>Bo 

Proof. By Proposition 12. 2[ it suffices to prove that for any differential operator L with 
constant coefficients we have 

(3.12) sup sup kB-^/^ sup I (LVb * 4) (2) I < 00. 

fc>0 B>Bo zeIR2 

Note that, by the standard symbol properties of (x)"'', we have 

\LVb{^)\<CVb{^), xGM^ 
where C depends only on Bq and L. Thus, it remains to prove that 

(3.13) sup sup kB-^^^ supKVb * 4)(^)| < 00. 

k>0 B>Bo zm^ 

We have 

VBiz) = {B-'\z\^ + l)-P/\ 



EIGENVALUE CLUSTERS 17 

Take z = {r,0), r > 0. Then 

1 r'^^ 

{VB*5k){z) = — / {B-\kcos9 - rf + B-\ksm9f + l)-p/^d9 
27r Jo 

2vr Jo Jo 



2 r'^ 

< 



TT 



/•TT/Z r) /"OO 

Jo Jo 

/ {{2e/7rY + iyp^^d0 = cB^/yk. 

Jo 



vrA; 

This yields 

sup \{VB*6k){z)\<CB'/^k~\ 
and f l333|) follows. □ 
3.3. Asymptotics of traces. 

Theorem 3.3. Let V G C^(M2). r/ien for each i eN, e>2, we have 
(3.14) Jim Af i)/2Tr(0p-(t^))' = |^ / { ^(c.,^^^^^^ 

The proof is based on the following technical lemma 



T JM. 



Lemma 3.4. Let i e N, i > 2, f e >S(M2{^-i))^ one? /ei the function ^ : x M^C^-i) ^ 
be given by 



(3.15) ip{iv,z) = ^Zj ■ (uj+i-Uj), 

where z = (zi, . . . , -2^_i) G M^*^^^"^-', uj = (ui, . . . , coe) G C M?^, and ■ denotes the scalar 
product in M^. Then 



(3.16) lim fc^-i / / /(z)e^^^('^'^)rfa;c/z 



= (27r)^ 1 / / f{aiu,a2Ui, ■ ■ ■ ,ai^iu) daida2- ■ ■ da^^idu. 
Jt Jk.^-^ 

Proof. The proof consists in an application of the stationary phase method. 
We use the following parametrisation of the variables a;, z: 

uji = (cos 6*, sin ^), 9 G [—71,71); 

Uj = {cos{6 + 0j), sm{6 + 6j)), 6j G [-vr, tt), j = 1, 1; 
Zj = ajUe + (ijOjf, ctj, /3j G M, j = 1, 1. 



18 



A. PUSHNITSKI, G. RAIKOV, AND C. VILLEGAS-BLAS 



We write « = («!,..., G R'~\ /3 = . . . , G = {6,, . . . , 6,,,) G 

[— TT, ttY~^. Using this notation, we can rewrite the integral in the l.h.s. of f l3.16p as 

(3.17) [ [ f{z)e"'^^'^'^Uu;dz 

F{a, 13, 6, 0)e^^'*("'^'^) df3 dct dO dO, 



-TV J (—77, 7r)^ 

where 

F{a, f3, 6, 9) = f{aiUJe + fiiu^, ae^iUi + A-icu/), 

and 

$(a,/3,0) = ai(l -cos^i) -/3isin^i, Hi = 2, 

$(a, /3, 6) = ai-i — {ai cos 6i + /3i sin 6i) 

i-i 

+ ^^((Q;j_i — aj) cos6j + {Pj-i — Pj) sinOj), if £ > 3. 

Let us consider the stationary points of the phase function $. By a direct calculation, 
V$(q:, 13,0) = if and only if /3 = and 6 = 0. By a standard localisation argument, 
it follows that the asymptotics of the integral f l3.17p will not change if we multiply F by 
a function x = x(/3, ^), X ^ C°°(R.^^^ x [— tt, tt)^"^), such that x(/3, 0) = 1 in an open 
neighbourhood of the origin f3 = 0, 6 = 0, and x(/3, 6) = if \f3\ > 1/2 or \0\ > 7i/2. 
Let us write 



(3.18) [ [ [ [ F{cx,f3,e,e)x{l3,e)e'^'^'^'^'^'^Ul3dcxdede 

I{k; a, 6) dct dO, 



where 



i{k-a,e)= f [ F{oL,j3,e,e)x{(3,e)e'^'^^'^^'^^^Uf3de. 



Let us fix OL, 9 and compute the asymptotics of the integral I{k; cx, 6) as k oo. A 
direct calculation shows that the stationary phase equations 

are simultaneously satisfied on the support of x if and only if /3 = 0, = 0. In order 
to apply the stationary phase method, we need to compute the determinant and the 
signature (i.e. the difference between the number of positive and negative eigenvalues) of 
the Hessian of $(«, /3, 0) with respect to the variables f3, 0. Let us denote this Hessian 
by H{a). The 2{i — 1) x 2{£ — 1) matrix H{ct) can be represented in a block form 



H{cx) 
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where 

An explicit computation shows that 

Hll = 0, H12 = {-Sg^p + ^Q.p+ljp^gLi , 

H21 = H22{a) = diag {ai, 0^2 - ai, • • • , a^-i - 0^-2} • 



Hence, 



det2(^_i)if(a) = det,_i(-ifi2i/2i) = (-1)'"'; 



in particular, our stationary point is non-degenerate. 

In order to calculate the signature of H{a.), note that since det H{a.) ^ for all cx and 
H[(x) depends smoothly (in fact, polynomially) on ex., the signature is independent of 
a. Some elementary analysis shows that sign H{0) = 0. 

Now we can apply a suitable version of the stationary phase method (see e.g. [T6| 
Chapter 1] or [121 Chapter III, Section 2]) to calculate the asymptotics of I{k](X,9). 
This yields 

lim k^-^I{k;cx,e) = (27r)^-^F(a, 0,0,0). 

Using, for example, the Lebesgue dominated convergence theorem, one concludes that 
the above asymptotics can be integrated over ex. and 6, see fl3.18p . This yields the 
required result f l3.16p . □ 

Proof of Theorem \3.3[ We use the notation t^, see (13. 4p . Denote by tk/ the Weyl symbol 
of the operator (Op'"(tfc))^, i.e. 

Op"(tM) = (Op"(^fc))^ ^ = 2,3,.... 

By the standard Weyl pseudo differential calculus (see e.g. [321 Chapter 7, Eq. (14.21)]), 
for Q G we have 

4/(0) = (27r)-^+i j 4(0 - 0-i)4(0-i - O-2) ■ ■ ■4(Ci)e^^^'=^"(^-^-^)(iC 
where cr(-, ■) is the symplectic form in x M^, and C, = (Ci, • • • , Ci-i)- It follows that 
(3.19) Tr(Op"'(t,))' = TrOp'"(tfe,^) = ^ / tfc,,(x)rfx = 4/(0) 

= (27r)-^+i / 4(-0-i)4(0-i-0-2)---4(Ci)exp(^ Va(0,0-i))^^C, 
where we use the convention that Y^j=2 = if £ = 2. 
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Recalling (jSl]), (jSZD, we get 

and so, substituting into fl3.19p . we get 

where 

/(C) = {27r)-'''-'VB{-Q-i)VB{0-i - 0.^2) ■ ■ ■ MCi) exp(i ^'ll ^iCv C-i)) , 
and !f is given by f l3.15p . Applying Lemma [3.4[ we obtain 

(3.20) \im k^-^Ti (Op'' {tk)Y 

= (27r)"^ / / VB{-ai^iUj)VB{{ai^i - 0^-2)^) ■ ■ ■ VB{{a2 - ai)uj)VB{aiUj) da du. 

It remains to transform the last identity into fl3.14p . We have 

VB{auj) = B ! e-'"''^'^'Vi{u,b)db 
Jm. 

where, in accordance with fl2.27p . we use the notation Vi{x, y) = V{—y, —x), (x, y) G M^. 
Therefore, 

(3.21) 

(27r)"M / VB{-ae-i(^)VB{{ae-i-ae_2)(^)---VB{{a2-ai)u)VB{aiu)da.duj 

5(i+^)/2 f f ^ 5(i+^)/2 r r ^ 

Viicu,bYdbduj = — / V{co,bYdbduj. 



Now (KWf follows from (K20\f and (ICTj) . □ 
3.4. Proof of Theorem 1.6. (i) First we observe that 

\\p,vp,\\ < iip,(r'^lii(r^iiii(r''/'nii 

< ii^iujin(r''/'iiii(r''/'^.ii = \\v\\xMP,{r'p,i 

and so it suffices to consider the case l^(x) = (x)~^. 
Next, by Corollary 12. 13[ we have 

\\P,VP,\\ = \\Op^{Vb*'^,)\\. 

In order to estimate the norm of Op"'(Vb * ^g), we use Lemmas 13.11 and 13. 2[ We note 
that for V{x.) = (x)~^, p > 1, we have V ^ and 

\viC)\<CN\C\-\ ICI>i, 
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for all > 1 (see e.g. [311 Chapter XII, Lemma 3.1]). Thus, the integral in the r.h.s. 
of (13. 2p is convergent, and so the proof of (13.21) apphes to ^(x) = (x)~^, p > 1. Now, 
combining Lemmas 13.11 and 13. 2^ we get 

sup sup AJ/25-1||0p'"(Vb*^,)|| 

g>0 B>Bo 

< sup sup Xf'B-'\\Op'"{VB * ^,) - Op'"{Vb * S^)\\ 

q>0 B>Bo 

+ sup sup X'J'B-'\\0p'"{Vb*6^)\\ < oo, 

q>0 B>Bo 

which proves the required estimate. 

(ii) As in the proof of part (i), we may assume V(x.) = (x)~^. First let us consider the 
case p > 2, i = 1. By Lemma [2.121 with £ = 1 we have 

B-'\\P,VP,h = ^ [ Vi^)d^<^\\V\\xJ (x)-''dx, 

which proves (II. lip in this case. 

Let us consider the case of a general i. For a fixed s > 1 and any i E [l,oo], let 

for £ = OO, one should replace l/£ by 0. By the previous step of the proof and part (i) 
of the theorem, 

sup sup II M^^) 111 <Ci<oo, sup sup ||M^°°)|| < Coo < oo, 

q>0 B>Bo q>0 B>Bq 

where the constants Ci, C^o depend only on Bq and s. Applying the Calderon-Lions 
interpolation theorem (see e.g. Theorem IX.20]), we get 

sup sup \\Mf%<cfct'^/'<oo 

q>0 B>Bo 

for all i > 1. It is easy to see that the last statement is equivalent to (II. lip . 

(iii) First let us note that the case £ = 1 is straightforward. Indeed, if the integer i = 1 
is admissible, i.e. if £ = 1 > l/(p - 1), then p > 2, V e L^R'^), and (12:291) yields the 
identity 

TrPgVPg = — [ V{x)dx = — [ [ V{Lo,b)dbduj. 
27r J^2 2-n Jj 

Thus, we may now assume i > 2. We will first prove the required identity (I1.12p for 
V G C^(]R^) and then use a limiting argument to extend it to all V G Xp. Denote 

(3.22) ^^(V) = ^[ [ V{uj,bYdbduj. 

Jj Jr 

By Corollary I2.13[ we have 

TiiPgVPgY = Tr(Op'"(VB * ^q)Y, 
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and Theorem 13.31 says 

lim Tr(Op'"(l^B * 5^))' = je{V). 

Thus, it suffices to prove that 

(3.23) hm AS^-i)/2|Tr(0p'"(VB * ^,)Y ~ Tr(Op'"(VB * '^vOT))'! = 0- 

q—>-co 

In order to prove (13.231) . let us ffist display an elementary estimate 

(3.24) iTr(A^i) -Tr(Ag| <£max{Pl||^M|A2||^^}||A-A2||,; 
here £ G N and An & Si, n = 1,2. The estimate follows from the formula 

A{-Ai = Y,4-'-\A,-A,)Ai 

3=0 

and the Holder type inequality for the 5'^ classes. 

Next, using Corollary 12.131 and part (ii) of the Theorem, we get 

(3.25) limsupA(^-^)/(2^)||0p'"(VB * ^g)||^ = limsup Af-^^/^^f) ||p^ypj|^ < ^_ 

q—^oo q—^oo 

Further, by estimate (13. 3p . using the assumption ^ > 2, we get 

(3.26) limsup A('^-i)/(2^)||0p"'(l^B * ^,) - Op^{Vb * 

< lim sup Aj/2||0p"'(rB * ^,) - Op'^iVB * 6^)h = 0. 

Combining ([321, and ^M), we obtain ^23\f for V e C^; thus, (11321) is proven 

for this class of potentials. 

It remains to extend (I1.12p to all potentials V E C(M^) that satisfy (II. ip . For £ > 
l/(p — 1), denote 

Ai{V) = lim sup A(^-i)/'Tr (P,VP,)^ 
6AV) = limm{Xi^-^^^^TT(P„VP„Y. 

q—>-oo ^ 

Above we have proven that 

(3.27) Ae{V) = 6e{V) = -fe{V) 

for all potentials V G C^(M^); now we need to extend this identity to all V G Xp. From 
(II. 5p we obtain, similarly to (I3.24p . 

|7^(V^i)-7^(V^2)| < 1^ / / \Vi{io,bY -V2{uj,bY\dbdu 

Jt Jr 
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It follows that 7^ is a continuous functional on Xp. Similarly, using f l3.24p and part (ii) 
of the Theorem, we get 

limsup Af-i)/2|Tr(P,ViP,)^ - Tr(P,V2P,)'| < Cmax{|| Vif^-/, || V2||^-;}||l^i - V^Wx,, 

q—^co 

and so the functionals A^, 6i are continuous on Xp. It follows that (I3.27P extends by 
continuity from to the closure Xp of in Xp. In order to prove f l3.27p for all 
V G Xp, one can argue as follows. For a given i > l/(p — 1), choose pi such that 
1 < pi < p and £ > l/(pi — 1). Then Xp C and by the same argument as above, 
dMZD holds true for all F G XO^. □ 

4. Proof of Proposition 1.1 and Theorem 1.3 
As already indicated, this section heavily uses the construction of [22j . 



4.1. Proof of Proposition 1.1. Set Ro{z) := {Hq — zl) ^. By the Birman-Schwinger 
principle, if A G ]R\U^Q{Aq} is an eigenvalue of the operator H , then —1 is an eigenvalue 
of the operator \V\^/'^RQ{\)V^f'^ . Hence, it suffices to show that for some C > and all 
sufficiently large g, we have 

(4.1) |||V^|V2i?o(A)|r|^/2|| < 1^ for all A G [Ag-P,Ag + S], |A-AJ>— . 



Choose m G N sufficiently large so that ||V^||/Am < 1/2, and write -Ro(A) as 

q+m ^ 

, '-^ Afc — A 

k=q—m 

Then, for A G [A, - B, Ag + 5], 



\mi^R,{X)\v\^i^\\< V "i*^ ^\mi^R,{X)\v\'i% 

, Afc — A 

k=q—m 

By the choice of m, one has 

\\\V\'/'Ro{X)\V\'/'\\ < \\\V\'/mi/XJ\\\V\'/'\\ = ||\/||/A™ < 1/2. 
On the other hand, by Theorem ll.6( i). 

lllT/lV2p |T/|1/2|| 

E ^'^ " <(2m + l)0(g-^/^) max \h - = Oiq-'/')\\ - \\-\ 

' ' A^j — A q—m<k<q+m 

k=q—m 

Thus, we get (14. ip for sufficiently large C > 0. □ 
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4.2. Resolvent estimates. Let Tg be a positively oriented circle of center Xg and radius 
B. 

Lemma 4.1. Let V satisfy (11.11) . Then for any i>l,i>l/{p — 1), one has 

(4.2) snp\\\V\^/^Ro{z)\V\^/^\U = 0{q-^'-^'>/^'\ogq), q ^ oo, 

(4.3) sup|||\/|^/^i?o(^)||2£ = Oiq-^'-'^/''\ogq), q ^ oo. 
Proof. Let us prove (14.21) . Using the estimate (11.111) . we get for z G Tg-. 

III I uv yi I ii« - ^ Xk-z \Xk-z\ 

- J, \B{2x + l)-z\ Jg^,\Bi2x + l)-z\ ^ 

= Oiq-^'-'y'nogq), 

as g — )■ oo. This proves (14. 2p . Using the fact that 

one proves the estimate (14. 3 p in the same way. □ 

4.3. Proof of Lemma 1.5. The fact that {PgVPgY e 5*1 follows directly from Theo- 
rem [L6l Let, as above, Tg be a positively oriented circle with the centre and radius 
B. Let q be sufficiently large so that (see Proposition II. ip the contour Tg does not 
intersect the spectrum of H. We will use the formula 

(4.4) (H - Xg)%x,-B,Xq+B)iH) = / (z- XgfR{z)dz, 

where R{z) = {H — zl)^^. Let us expand the resolvent R{z) in the r.h.s. of (14.41) in the 
standard perturbation series: 

oo 

(4.5) R{z) = Ro{z) + J2i-^yMz){VRo{z)y. 

Let us discuss the convergence of these series for z G F^, g large. Denote W = IV]^^"^, 
Wq = sign(y). For j > £, we have 

(4.6) \\Ro{z){VRo{z)yh = URo{z)W){WoWRo{z)Wy~'Wo{WRo{zm, 

< \\Ro{z)W\\2i\\WRo{z)W\\r'\\WRo{z)he, j > t 

Applying Lemma 14. H we get that the series in the r.h.s. of (14.50 converges in the 
trace norm for z & Tg and q sufficiently large (note that although the tail of the series 
converges in the trace class, the series itself is not necessarily trace class). 
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Next, it is easy to see that the integrals 

iz-XgYMz){vMz)ydz 

with j < i vanish (since the integrand is analytic inside Tg). Thus, recalling f l4.6p . we 
obtain that the operator {H — \qYl(^x^^B,\q+B){H) belongs to the trace class and 



(4.7) Tr{(i7-A,ri(,, 



1 °° r 

j=e -^^1 



z)y]dz. 



Integrating by parts in each term of this series and computing the term with j = i hj 
the residue theorem, we obtain 

(4.8) Tr{(i7-A,)^1(,,_B,A,+B)(i/)} 

= Tr(P,VP,)^ + ^ E ^ / {z-X,Y-'MVRo{zWdz. 
j=i+i -^^1 

It remains to estimate the series in the r.h.s. of f l4.8p . This can be easily done by using 
Lemma [4.11 Similarly to (14. 6p . we have 

\TT{VR^{zy^ < \\WRo{z)W^ < \\WR^{z)W\\l, j > i, 



and so 



[z- \gY~^TT{VRo{zy^dz 



<C,iC,q-^'-'y''\ogqy, 



for all sufficiently large q. Thus, the series in the r.h.s. of (14.80 can be estimated by 

oo 

Ci Yl Ciq-'^^ilogqy. 

j=i+i 

For all sufficiently large q, this series converges and can be estimated as o{q~^^~^^^'^). □ 

4.4. Proof of Theorem 1.3. Let R > Ci where Ci is the constant from Proposi- 
tion [TTTl Then 

(4.9) ll[_«,«](Af (i/-A,)) = ll(,,_B,A,+B)(i^), qeZ+. 

Next, choose -R > Ci so large that suppf? C [—R, R]. Let Eq be an even natural number 
satisfying > l/{p — 1). Since g{X) by assumption vanishes near A = 0, the function 
q{\)/\^° is smooth. Applying the Weierstrass approximation theorem to this function 
on the interval [— i?, i?], we obtain that for any e > there exist polynomials P+, P_ 
such that 



(4.10) 
(4.11) 



P±(0)=P4(0) 



P 



(4-1) 



(0) = 0, 



P_(A) < ^(A) < P+(A), VAg[-P,P], 
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(4.12) P+ ( A) - P_ ( A) < eX^" , VA G [-R, R] . 
Thus, we can write 

for any A G [— -R, R], and therefore 

Tr{l[_^,«](Af - A,)) P-{K^'{H - A,))} < Tr ^.(Af - A,)) 

By (14. 9 p it follows that for all sufficiently large q, 

(4.13) Tr{l(,,.^,,^+^)(if) P-(Af (i^ - A,))} < Tr f,(Af (iJ - A,)) 

< Tr{ll(,,_B,,,+^)(i7) P+(Af (i^ - A,))}. 
By Lemma [1.51 and Theorem ll.6( iii). we have 

hm A-^/2 TT{l(^x,-B,x,+B)iH) P±{\\'\H - A,))} 

q^oo 



1 

2^ 




P±{V{u,h))dhduj = I P±{t)dfi{t). 

Jr 



Combining this with (I4.13p . we get 

limsupA^i/2Tr^(Aj/2(i/- Aq)) < [ P+(A)d/x(A), 

limM X-^/^Ti g{Xl/\H - X.)) > [ P_(A)ci/i(A). 

Finally, by (KWf . 

'{P+iX) - P_(A))rf/i(A) <e [ A^«rf/i(A). 



By (II. 6p . the integral in the r.h.s. is finite. Since e > can be taken arbitrary small, 
we obtain the required statement. □ 

Appendix A. 
A.l. Proof of formula (3.6). By definition, 

(see dsn). Further, by [H Eq. 22.12.6] we have 



m=0 
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where Lm^^'^\ m G Z+, are the generahzed Laguerre polynomials of order —1/2. By [H 
Eq. 22.5.38] we have 

m!/ 

where are the Hermite polynomials. Therefore, 

T {\ |2n _ H2m(Ml)H2g-2m(^2) 

HU^I J - 2^ m!(g-m)! ' 

m=0 ^ 

and 

2^^+^7r^ ^^m\{q - my. 

It is well known that the functions Hm(t)e~*^/^, t G M, m G Z^., are eigenfunctions of 
the unitary Fourier transform with eigenvalues equal to (— i)"^ (see e.g. [5J). Hence, 

^(0 = ^ E u ^ „ H2^(2-^/^Ci)H2,-2^(2"^/^C2)e-l^l^/^ = 

m=0 ^ 

(27r)-iL,(2-i|Cne-lfl'/^ = (-l)^vI/,(2-iC)/2. 

A.2. Proof of estimate (3.10). Denote Uq{x) = e'^'/^Lqix), Vq{x) = Jo( V(4g + 2)x). 
Using the differential equations for the Laguerre polynomials and for the Bessel func- 
tions, one easily checks that Uq and Vq satisfy 

XUg{x) + U'q{x) + (g + \)Uq{x) = jUq{x) , 

xvg{x) + v'g{x) + (g + ^)vq{x) = 0. 

Using these differential equations and the initial conditions for Uq{x), Vq{x) at x = 0, it 
is easy to verify that Uq satisfies the integral equation 

(A.l) Ug = Vq + KqUq, {Kqf){x)= Fg{x , y) f {y)dy , 

Jo 

Fq{x, y) = -jy{U^m + 2)x)YoWm + 2)y) - Y,{^ + 2)x)JoW {Aq + 2)y). 
This argument is borrowed from [SU]. Iterating f lA.ip . we obtain 

(A. 2) Uq-Vq = KqVq + K^Uq. 

Now it remains to estimate the two terms in the r.h.s. of flA.2p in an appropriate way. 
Using the estimates 

|Jo(x)| <C/v^, |lo(x)| <C/v^, x>0, 

we obtain 

\Fq{x,y)\<Cq-^''^x-^'S'^'^. 9 e N, z > 0. 
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This yields 

(A.3) r F^{x, y)v,{y)dy < Cq-^/'x-'" f y^'^dy = Cq-'/'x'>'\ 

Jo Jo 

Next, using the estimate < 1 (see Eq. 22.14.12]), we obtain 

\iK,u,)ix)\ < Cq-"^x-^'^ r y^'^dy = Cq-"^x^l\ 

Jo 

and so 

(A.4) \{Klu,){x)\ < Cq-^x-^l^ Py^^-^dy = Cq-^x\ 

Jo 

Combining (1A.2P with (1A.3P and (IA.4p . we obtain the required estimate ( I3.10p . 
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